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I. INTRODUCTION 



In a previous article [1], we have considered a simple three-flavor neutrino-oscillation 
model with both mass-square differences (Am--) and timelike Fermi-point splittings (Ab^). 
The mixing of the mass sector is taken to be bi-maximal and the one of the Fermi-point- 
splitting sector trimaximal, with all complex phases vanishing. The model has furthermore a 
hierarchy of Fermi-point splittings (b^ = b^ ^ b^) which parallels the hierarchy of masses 
(mi = m 2 7^ m 3 ). This particular model may be called the "stealth model," as it allows for 
Lorentz-noninvariant Fermi-point-splitting effects to hide behind mass- difference effects. For 
the physics motivation of this type of model discussion, see Ref. [1] and references therein. 

The present article follows up the previous one in three ways. First, we give model results 
of the vacuum appearance probability P{y^ — > u e ) relevant to MINOS in the low-energy 
mode, complementing the specific results for the medium energy mode in Ref. [1]. 

Second, we consider the model probability P^ e = Piy^ — > u e ) for the case of relatively 
strong Fermi-point-splitting effects compared to mass-difference effects, whereas Ref. [1] 
focused on relatively weak splitting effects. In particular, we introduce a new parametriza- 
tion (with nonnegative dimensionless parameters p and r) which makes a straightforward 
comparison between different long-baseline neutrino-oscillation experiments possible. Rel- 
atively weak or strong Fermi-point-splitting effects then correspond to pr <C 1 or pr ^ 1, 
respectively. The behavior of P^ e {p, r ) turns out to be quite complicated for pr ^ 1, 

Third, we present results on the appearance probability Piy^ — > v e ) from a generalized 
model with the same mass hierarchy as the model of Ref. [1] but with equidistant Fermi- 

(2) (1) f3) (2) 

point splittings (% — % — % % ) and one additional complex phase (cu = 7T/4). For 
completeness, we also give the model probability of the time-reversed process, u e — > v^. 
It will be seen that the generalized model has a rather interesting phenomenology with 
stealth-like characteristics in certain cases and strong time-reversal noninvariance in others. 

The outline for the remainder of this article is as follows. In Sec. II, we recall the definition 
of the model and introduce the two neutrino-oscillation parameters p and r. In Sec. Ill, we 
present and discuss our new results. In the Appendix, we consider one concrete model and 
compare the capabilities of four accelerator-based long-baseline oscillation experiments. 

II. MODEL 

Setting H — c — 1 and writing p = \p\ for the neutrino momentum, the Hamiltonian of 
the stealth model [1] contains three terms in the (v e , v T ) flavor basis, 

HD P l + X-D m -X- 1 +Y-D bo -Y-\ (1) 
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with diagonal matrices 



D m = diag ( g , g , g ) , A.« diag ( 6<" , tf> , tf> ) 



(2) 



and S77(3) matrices 

X = M 32 (# 32 ) • M 13 (0 13 , 5) ■ M 21 (0 21 ) , ^ = M 32 (x3 2 ) • M 13 (xi 3 , • M 21 ( X2 i) , (3) 

in terms of the basic matrices 

I I \ / costf sintf o\ 

M 32 (tf) = costf sintf , M 21 {&) = 

^ — sin?? cos-# y 

+e+^sintf\ 



— sin i? cos 
\ 1 



M 13 (0,<p) = 



and the following parameters: 



^ cosfi 

1 

^-e-^sintf cos$ y 



A< =ml-ml = 0, R = Ab^/Ab^ = (tf> - b ^)/(bf - = 



e 



13 



0, : 



21 







32 



Xl3 = X21 = X32 = 7T/4, 



(4) 

(5a) 
(5b) 



<5 = = . (5c) 
With all other complex phases vanishing, there are only two neutrino-oscillation parameters 
left in the model, 



Am 2 31 



m 3 - m( > , Ab ( 



(31) 



l(3) 



^>o, 



(6) 



which have been taken positive. 

For high-energy neutrino oscillations over a travel distance L, there are, then, two di- 
mensionless parameters which completely define the problem, at least for the simple model 
considered and matter effects neglected. These neutrino-oscillation parameters can be de- 
fined as follows [E v ~ p): 

2E„hc f 2.5 x 1(T 3 eV 2 /c 4 ^ ^735 km\ 



P = 



(Am^l c 4 L 



\Am 2 31 \ 



L J \ 4.656 GeV 



L\Ab, 



(31), 



IA6, 



(31) | 



L 



(7a) 



(7b) 



he V 2.685 x 10- 13 eV I \735 km, 

with h and c temporarily reinstated. 

In terms of these parameters, an approximate formula for the vacuum probability P^ e = 
P(v^ -> v e ) is given by [1] 



prteaith (p>r) _ (1/2) sin 2( 2 13 ) sin 2 Q p -l + r + ^-2 + r 2 ] /4) , 

with the following energy-dependent effective mixing angle: 
9i 3 (/9,t) ~ (1/2) arctan(pr) . 



(8a) 



(8b) 
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Recall that, according to Eq. (5b), the bare mixing angle 9 i3 vanishes identically and that, 
as mentioned above, matter effects are neglected in this approximate model probability [2]. 
More generally and allowing for a nonzero Fermi-point-splitting ratio R, one can define 

P; e tcalth (P, t) = sin 2 (© 32 ) sin 2 (2 13 ) A 2 e , (9a) 

with the following functional dependence: 

e 32 = §32 (pr,R), 13 = e 13 (pr, J R), V = A Me (p , r , R) e [-1, +1] , (9b) 

so that the effective mixing angles O32 and 613 are independent of the travel distance L, 
according to Eqs. (7ab), whereas A^ e does depend on L [5]. In principle, 613 does not have 
to vanish in the limit p — > for fixed r and this is indeed the case for the generalized model 
with splitting ratio R 7^ (see below). 

III. RESULTS AND DISCUSSION 

In Ref. [1], we have given five figures with model results, three of which may be relevant 
to T2K or NOz/A and two to MINOS in the medium-energy (ME) mode. Here, we present 
one more figure, Fig. 1, which may be relevant to MINOS (baseline L = 735 km) in the 
low-energy (LE) mode, with peak neutrino energy E u approximately equal to 3.75 GeV and 
neutrino-oscillation length scale L = 2ixE u /\Am\ l \ ps 1860 km for Am^ = 2.5 x 1CT 3 eV 2 . 
Figure 1 (which may be compared to Figs. 4 and 5 of Ref. [1]) shows that, if MINOS-LE 
would be able to place an upper limit of 10 % on the appearance probability Piv^ — > v e ) 
at the high end of the energy spectrum [E v > 4 GeV), this would correspond to an upper 
limit on |A&q 31 ^| of approximately 3 x 10~ 13 eV (assuming r < 7r, see below). 

A further figure, Fig. 2, gives numerical results which illustrate the general behavior of the 
vacuum probability P^ e = Piy^ — > v e ) as a function of the two dimensionless parameters p 
and r from Eqs. (7ab). Note that we expect lim r ^ P^e(p, t) to vanish [pure mass-difference 
model with Am^ = and 9 13 = 0] and lim^oo P^ e (p, t) to be given by (1/2) sin 2 (r/2) 
[pure Fermi-point-splitting model with Ab^ = 0, Ab^ 7^ 0, and trimaximal mixing]. The 
landscape of Fig. 2 can then be described as follows: mountain ridges start out at p 3> 1 and 
r ~ rioo Ti for odd integers noo, slope down towards lower values of p keeping approximately 
the same values of r, and, finally, disappear while bending towards lower values of r (more 
so for ridges with large n^). This topography is qualitatively reproduced by the analytic 
expression (8). 

Figure 3 presents a sequence of constant-r slices of the vacuum probability P^eip, t) from 
Fig. 2. The behavior of P Me (p, r) at r = 2tt (or integer multiples thereof) is quite remarkable, 
being nonzero only for a relatively small range of energies; cf. the short-dashed curve in the 
upper right panel of Fig. 3. 

Next, turn to a generalized stealth model with splitting ratio R = 1 (i.e., equidistant 
Fermi-point splittings) and complex phase u = or it/ 4, the other model parameters being 
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kept at the values (5abc) and maintaining the signs (6). (For high energies, this model is 
similar to the pure Fermi-point-splitting model studied previously [6, 7].) Figures 4 and 5 
give the resulting vacuum probabilities P^ e = Piy^ — > is e ) [8, 9]. For both complex phases, 
the behavior of P Me (p, r) at r m 12 is particularly noteworthy; cf. the short-dashed curves 
in the lower right panels of Figs. 4 and 5. At the corresponding distance L for given value of 
the Fermi-point splitting A6q 31 \ namely, the stealth model lives up to its name by evading 
detection via v e appearance, unless the experiment is able to reach down to low enough 
neutrino energies (p ~ 0.15). In principle, the way to corner this stealth model would be to 
use several broad-band experiments at different baselines, but this may require a substantial 
effort (see the Appendix for a case study). 

In the previous paragraph and corresponding Appendix, the stealth-like behavior of the 
R — 1, uj — 7r/4 model has been emphasized, but this holds only for the z/ M — > v e channel 
relevant to superbeam experiments with an initial beam from pion and kaon decays. In 
other channels, the situation may be different. The probability P efl of the u e — > channel, 
for example, is given in Fig. 6. The very different probabilities of Figs. 5 and 6, for p > 0.2 
and generic values of r, signal time-reversal (T) noninvariance. 

For comparison, we give in Figs. 7-10 the numerical results from a "realistic" mass sector: 
R m = Am^/Am^ ~ 0.033, sin 2 (2#i 3 ) = 0.1 [close to the experimental bound from Chooz], 
and two possible values of the complex phase, 5 = or 5 = ir/2. Figures 9 and 10, in 
particular, show that the maximum T-violation discriminant A T = P efJl — P Me from pure 
mass-difference neutrino oscillations (r = 0) is of the order of several percent, much less 
than the potential Fermi-point-splitting result (several tens of percent) at the high-energy 
end of the neutrino spectrum. 

The strong high-energy T violation (and possibly CP violation [7]) from Figs. 5-10 traces 
back to the large complex phase u of the model considered, together with the large mixing 
angles Xij an d splitting ratio R in the Fermi-point-splitting sector. In other words, the 
breaking of time-reversal invariance would primarily take place outside the mass sector. A 
neutrino factory (see Ref. [3] and the Appendix) would be the ideal machine, in principle, 
to establish such strong T violation in high-energy neutrino oscillations. 

ACKNOWLEDGMENTS 

It is a pleasure to thank Jacob Schneps for useful discussions and Elisabeth Kant and 
Christian Kaufhold for help with the references and the figures. 

APPENDIX A: CASE STUDY 

In this appendix, we give a concrete example of the complicated phenomenology of 
stealth-type models. Take, then, the generalized model (l)-(6) with parameters uj = 7r/4 
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and A&0 31 " 1 = 2A&Q 21 - 1 = 3.5 x 1CT 12 eV and consider the combined performance of four 
accelerator-based neutrino-oscillation experiments [10, 11]: the current K2K experiment 
(baseline L = 250 km and peak energy E u ~ 1.0 GeV), the running MINOS experi- 
ment (L = 735 km and E v pa 3.75 GeV, in the LE mode), the planned T2K experiment 
[L = 295 km and E v pa 0.7 GeV, for an off-axis beam at 2 degrees), and the proposed NOz/A 
experiment (L — 810 km and E v pa 2 GeV, at 14 mrad offset in the ME mode). At the end 
of this appendix, we also comment briefly on the capabilities of a possible neutrino factory. 

First, calculate the model predictions for the two current experiments. For K2K, the 
neutrino-oscillation parameters would be r pa 4.4 and p pa 0.63, so that a value for 
P Me = P{y^ — > v e ) of only 1 % would be expected (cf. solid curve in upper right panel 
of Fig. 5), which is consistent with the experimental result [13]. [Note that the K2K appear- 
ance probability P Me would be approximately 15 % for the model with vanishing complex 
phase uj, as indicated by the solid curve in the upper right panel of Fig. 4.] For the MINOS- 
LE experiment with r pa 13.0 and p pa 0.81, a relatively small probability P^ e pa 7% would 
be expected, which would perhaps be hard to separate from the background. [Note that this 
r-value corresponding to the MINOS baseline is an order of magnitude above the largest 
one of Fig. 1.] Hence, the model predictions of the appearance probability P Me , for the 
chosen parameters, would be consistent with rather low v e event rates from both of these 
current experiments. For completeness, the model would also give a survival probability 
P w = -> of some 49 % for K2K and 87% for MINOS-LE. 

Next, turn to the model predictions for the two future experiments considered. For T2K 
with r pa 5.2 and p pa 0.38, a substantial P^ e of some 25 % would be expected (cf. long-dashed 
curve in upper right panel of Fig. 5). Similarly, for NOz/A with r pa 14.4 and p pa 0.39, P^ e 
would be approximately 27%. Hence, the model predictions of the appearance probability 
P Me would imply a clear signal in these future experiments, at least for the chosen parameters. 
Again, for completeness, the model survival probability P MM would be approximately 71 % 
for T2K and 20 % for NOz/A. 

Ultimately, the lowest values of IA&q 31 - 1 ! could be probed by a neutrino factory with broad 
energy spectrum E v pa 10 — 50 GeV and several detectors at baselines L up to 12800 km; 
see, e.g., Refs. [3, 4, 12]. For a quick estimate, one can simply compare the expected 
experimental capabilities with the approximate vacuum probability (8), for p = 0(1) and 
r < 1 as defined by (7). In the best of all worlds, a P Me sensitivity of 0.01 % at L = 7350 km 
and E v pa 46 GeV, for example, would correspond to lAfeg 31 ^ of the order of 10~ 15 eV. 
However, for a definite analysis at these large distances, matter effects [2] would need to be 
folded in. 
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FIG. 1: Left: Fixed-distance vacuum probabilities P^e = P{y^ — > u e ) at I = L/L = 
L\Am^ 1 \/(2TvE l/ ) = 735/1860 vs. dimensionless energy x = E v jE v for the "stealth" model of 
Ref. [1] with both mass differences and timelike Fermi-point splittings, here defined by Eqs. (1)- 
(6). The broken curves are given by the analytic expression (4.4) of Ref. [1] with parame- 
ters Gig = ~E u \Abf l) \/\Aml 1 \ = 0.15, 0.30, 0.45, which correspond to Fermi-point splittings 
\Ab { Q 3) \ w (1.0, 2.0, 3.0) x 10~ 13 eV for E v = 3.75 GeV and Am^ = 2.5 x 10" 3 eV 2 . The 
same analytic expression is given here by Eq. (8) with parameters r ~ (0.37, 0.74, 1.12). These 
broken curves are only approximate and the corresponding thin solid curves give the full numerical 
results. The heavy solid curve gives, for comparison, the standard mass-difference probability -P^ D 
for Am|i = 0, sin 2 ^32 = 1/2, sin 2 (26>i 3 ) = 0.2, and I = 735/1860. Right: Electron-type neutrino 
energy spectra f e {x) at I = 735/1860 from an initial muon-type spectrum (4.8) of Ref. [1] at I = 
and vacuum probabilities P{v^ — > v e ) of the left panel. The initial spectrum multiplied by a 
constant factor 0.25 is shown as the heavy long-dashed curve. The resulting v e spectra from the 
"stealth" model are shown as the thin solid curves [approximate values as thin broken curves]. 
The heavy solid curve gives, for comparison, the v e spectrum from the standard mass-difference 
probability P^jP of the left panel. 
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FIG. 2: Numerical results for the vacuum probability P^ e = P(fu — > v e) from the "stealth" model 
(l)-(6), as a function of the two dimensionless parameters p and r, defined by Eqs. (7ab). Left: 
surface plot. Right: contour plot, with equidistant contours at P^ e = 0.05,0.10,0.15, . . . 0.45 and 
contours at integer multiples of 0.10 shown dashed. 
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FIG. 3: Numerical results for constant-r slices of the vacuum probability P = P(yu - > v e) from 
Fig. 2. The curves for positive r = 1,2,0 (mod 3) are shown as solid, long-dashed, and short- 
dashed lines, respectively. For the simple model considered, given by Eqs. (l)-(6), the parameters 
R and uj vanish identically. 
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FIG. 4: Same as Fig. 3 but for the generalized model with splitting ratio R = |A6q /A6q | = 1, 
the complex phase uj still vanishing. 
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FIG. 5: Same as Fig. 4 but now for the generalized model with complex phase u = vr 
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FIG. 6: Same as Fig. 5 but for the time-reversed process, with probability P' = P(v e 
CPT invariance holds, P' also corresponds to P(V^ — ► V e ). 
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{R m , sin 2 (20 13 ), 5 ; R, sm 2 (2 Xl3 ), oj} = {1/30, 1/10, ; 1, 1, tt/4} 
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FIG. 7: Same as Fig. 5 but now for mass-sector parameters R m = Am^i/Am^ = 1/30, 
sin 2 (2^i3) = 1/10, and (5 = 0. The heavy solid curve in the upper left panel has r = 0, cor- 
responding to pure mass-difference neutrino oscillations. 
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FIG. 8: Same as Fig. 7 but for the time-reversed process, with probability P' = P{v e 
CPT invariance holds, P' also corresponds to P(pu — ► V e ). 
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FIG. 10: Same as Fig. 9 but for the time-reversed process, with probability P' = P(v e 
CPT invariance holds, P' also corresponds to PiV^ — > z7 e ). 
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